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Abstract 



We get back to the computation of the leading finite size corrections to some random link 
matching problems, first adressed by Mezard and Parisi [J. Physique 48 (1987) 1451-1459]. 
In the so-called bipartite case, their result is in contradiction with subsequent works. We 
show that they made some mistakes, and correcting them, we get the expected result. In 
the non bipartite case, we agree with their result but push the analytical treatment further. 

1 Introduction 

The possibility of investigating some optimization problems using techniques of the field of 
disordered systems in statistical physics has been recognized for a long time (see p[ for a recent 
review). Along with the traveling salesman problem and, more recently, K-Sat, one of the 
problems which got most of the attention is the matching problem, of which two variants have 
been studied: 

(a) The simple matching problem: consider 2N points and a set of 'distances' between them 
lij = Iji. A matching of these points is a set of N pairs so that each point belong to exactly 
one pair. The 'length' of such a matching is 



One focuses on the properties of the matching of minimal length. 

(b) The bipartite matching problem (or assignment problem), which is as above, except that 
we split the points into two sub-sets A and B of N points each, and allow only matchings 
where each pair is made of a point from A and a point from B. 

Here we are interested in the case where the iy are independent identically distributed ran- 
dom variables, either uniformly distributed on the interval [0, 1] (so-called Eat case) or distributed 
with the law exp(— I) on [0, +oo[ (so-called exponential case). 

Both the simple and bipartite cases have been investigated in the thermodynamical limit 
N — * +oo, where self-averaging of the optimal length occurs. The replica method, in the replica 
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symmetric scheme, yielded predictions for the mean optimal length and the distribution of the 
lengths of occupied links in the optimal configuration |J. This was shown to be equivalent to 
a cavity approach Q . Numerical works checked the validity of the results obtained with these 
techniques ||] , and got interested in another quantity, the probability for some given point to 
be connected to its A;-th nearest neighbor in the optimal matching. This was also dealt with by 
an analytical cavity computation Q. Remarkably, |5J confirmed all the above results by rigorous 
proof. 

The stability of the replica symmetric solution was checked in Q, yielding as a byproduct 
the 0(1/N) correction to the mean length of the minimal matching. For the assignment problem 
with flat distances, they found 

i£ = £-i(g + 2 C (3)) +0 (I), (2) 

where 7TT means the average with respect to the distribution of the distances. 

This seemed to agree with the numerical simulations at the time Q . But || came up with a 
conjecture for the assignment problem at any finite ./V: in the exponential case the mean length 
of the optimal matching would be 

N , 

jexp I _ J_ (o\ 
^min\N / t ^ 2 ' V ' 

k=l 

which implies 

LZL = \-^ + o(^). (4) 
The problem is that in the framework of JU it is not difficult to prove - as we shall - that 

-jexp _yjlat _ 2C(3) 

so that (HJ) is not compatible with (||). Beside, other results sustain (0): more recent numerical 
simulations , and an allegedly rigorous proof of (||) . 

The main purpose of this article is to show that Mezard and Parisi [|| actually made some 
mistakes in the computation leading to @. Correcting them, one does get 

^ = y-^(l + 2C(3))+o(l), (6) 

which, using (0), coincides with (||). 

Apart from some trivial errors due to a confusion about the number of points, iV or 2N, 
that has already been pointed out Q, and some typos, there are essentially two mistakes in ||]: 

(i) they forgot a contribution to the 0(l/iV) term of L m i n , but this is without any consequence 
as this term turns out to vanish in the zero temperature limit (appendix |A.1| ); 

(ii) they made a mistake in the computation of another contribution (appendix |B.3j ), and this 
was responsible for the wrong result 

As the computations carried out in Q are quite involved, we have chosen to make this article 
the most self-contained possible by restating all the necessary steps. 

In section 0, we tackle the simple matching problem, which is formally simpler but very 
similar to the assignment problem. Error (i) is common to both problems, and we deal with it 
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in that section. Moreover we refine the computation of the 0(1/N) correction which in [|| relies 
on a rough numerical procedure. We give an analytical expression of the correction as the sum 
of a series. Unfortunately we were not able to sum this series, but it might not be impossible. 

In section ^, we turn to the assignment problem. There we correct error (ii), which is specific 
to this variant. 



2 Non bipartite case 

In this section, for the sake of simplicity, we exclusively consider the Bat case. 

To tackle the problem with the tools of statistical mechanics, one introduces an inverse 
temperature f3 (to be sent to +00 in the end) and a partition function 

Z = E exp(-N(3L matching ). (7) 

all possible matchings 

The scaling factor N ensures a good thermodynamic limit at fixed /3 [|l^] . 

We will not get into the details of the computation of the averaged replicated partition 
function, because it is quite similar to the bipartite case, for which the derivation can be found 
in appendix |B.1|. Let us just state the result H: 



^=/n^=-p(-f^[Q])expf-itff I ^u 7 V (8) 

J „ . Inn-.,.* N V 2 J \ 4^ff aUl '/ 



where g p = l/(f3p) and 

S [Q] = J2 — -^*[Q], (9) 



g. 



*[Q] 



In the above expressions, a stands for any non empty subset of {1, ... n} so that the number 
of Q a variables is 2 n — 1. For such an a, we call p(a) its cardinal number, and use the shorthand 
notation g a for g p ( a y The notation Y^'ocy m eans that the summation runs over all the couples 
(a, 7) such that a n 7 = 0. 

Lmin is nothing but the intensive free energy F/N = —l/((3N)h\Z in the limit (3 — * +00. 
It is evaluated by a saddle-point method. The saddle-point equation reads 

^=4^. do) 

g a dQ a 

It has been solved under the assumption of replica symmetry Q: Q S J> = Q S p P ( a y It turns out 
that the order parameters are not well defined quantities at low temperature, and one can 
bypass this difficulty by considering the well defined generating function 



pi 

p=i 1 



for which ( JlOj ) translates into 

au\ = - 

13 



G(0 = -| / dyK{l + y)e- G{ -y\ (12) 
00 
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+°° (— 1 \p 

K(u) = J2 -^T ePU = - 1 + Jo(2e«/ 2 ). (13) 



where 



In ( |13|) Jo is the Bessel function of order zero. Note that (|12| ) can be obtained by direct 
probabilistic arguments using the cavity method [E|. 

The free energy in the thermodynamical limit has been evaluated as 

Jf = & 2]M SlQ " ] = (14) 

/3^+oo 

The first finite size correction AF is the sum of two terms: the first one, corresponding to 
the last factor in (||), 

has been computed in the appendix of [H. The second one corresponds to the Gaussian fluctu- 
ations around the saddle-point: 



AF 2 = — lndet 

2f3Nn 



1 d 2 S 



2 V9«V9, dQaQ 



' lndetT+ 7 , (16) 



2/3 Nn 

where we define the matrices T^ 7 , e = dbl, as follows: 

nspfisp qsp 

rre x , ^« V7 x ^aU7 , /-, 

i a7 — d a~f + £ o . = ~ g^an 7 =0 \/9a9j- {<-<) 

At this stage the introduction of e is a useless complication, since for the non bipartite case 
we only have to consider the case e = +1. However the case e = — 1 will turn out to be useful 
in section |3], where we deal with the bipartite case. 

The computation of AF 2 involves finding the eigenvalues of (T^ 7 ). A vector (Q a ) is eigen- 
vector for the eigenvalue A if 

Q sp Q sp Q sp 
Va, Q a + £ Y] Qi T T, 7 ~ £ Q 7 -v/gaffr = ( 18 ) 

7 V&*97 an7=0 9aUy 



The diagonalization process adopted by [g[ follows the de Almeida-Thouless strategy ]13|] of 
considering stable subspaces of increasing dimension. First we look for eigenvectors where no 
particular replica is distinguished: 

Qa = C p{a y (19) 

When we plug this into (|l8|), we see that we have to diagonalize a n x n matrix N°' £ (n) 



n s P nspnsp 

^P+q r—— , yh 

J n.-p _ Vi/py? o /7r - 7r 

5p+<? z \J9p9q 



N°> e (n) pq = S pq - sCl_ p ^Jg-g- q + e^^C* p,q = l,...n, (20) 
the 

N°' £ when n — > 0: 



where we use the notation C£ = ,/"' w . This matrix turns into an infinite dimensional matrix 
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, (p + q-l)\ Q SP 



< £ = S pq - ^ p+q Va^ P,g = l,... + oo (21) 



(the last term in (BO) does not contribute). 



The eigenvalues of N ' 6 have multiplicity 1 in the spectrum of T £ 

- «7 



Then we look for eigenvectors of T £ where one replica is distinguished, say a: 



Qc 



dp( a ) if a G a 
e p ( a ) if not 



The orthogonality constraint between this family and the previous one (119) reads pd p + (n — 
p)e p = 0. So we can choose the only d p as variables (of which d n = 0), and we end up with the 
diagonalization of an (n — 1) x (n — 1) matrix N 1,£ (n) 



(l sp n n s Pn s P n n s Pn s P 

p,g = l 5 ...n- 1. 



We get n—l eigenvalues Ai, . . . X n -i (independent on a) and the corresponding eigenvectors 
uf, . . . u®_ l . The important point is that the uf, a = l,...n are not linearly independent: 
Yl a u i ^ s a vec t° r °f the previous family and orthogonal to it, so it has to be 0. Eventually the 
eigenvalues of N 1,e (n) have multiplicity n — 1 in the spectrum of T^ 7 . When n — > 0, iV 1,e (n) 
becomes an infinite dimensional matrix which happens to be exactly N°' £ . 

More generally, one finds the whole spectrum of T^ 7 by looking for eigenvectors which have 
k given distinguished replicas: 



Qa 



if p < k 

dp(a) ^ a con t a i ns k + 1 — i of the distinguished replicas, 



where i goes from one to k + 1. The orthogonalization with respect to a family where only k — 1 
of these replicas are distinguished reads as a system of equations 

Vi = 0, 1, . . .k - 1, J2 Cl_fl-_t +j) d k P +1 ~ {r+3) = 0> (23) 

r=0 

whose solution in the n — > limit is 



4 4 



p(p + 1) . . . (p + fc - 1) (p - A; + l)(p + 1) . . . (p + k - 1) 



(24) 



<ip dp - * - 



(p-A; + l)(p-fc + 2)(p + 2)...(p + fc-l) (p-A; + l) 



(note that this is slightly different from equation (20) in || where there is a typo) 

It follows that we can keep the only dl as independent variables, and have to diagonalize a 
matrix N k,£ (n), which in the limit n — ► is the infinite dimensional matrix 

- s p „ - 4-y fr 1V f^VW, M -i,... + =. (25) 
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The eigenvalues of N k,£ (n) have multiplicity C k — C k 1 in the spectrum of . 

The article || prefers using matrices derived from the N k,£ by some transformations which do 
not affect the spectrum: after shifting the indices p and q by k, then transposing and multiplying 
each entry by (— l) p+q * I , one gets the family of matrices 



9 P +k (p+l)...(p+fc-l) > 

*t = L - g (-l)^ (P + g + 2fc " 1)! ^^^. 



P ' q {p + 2k-l)\q\ g P + q+ 2k 



where p, q = 0, 1, . . . + oo. 

We can now proceed to the computation of AF 2 ((If 



AF 2 = AF 2 ' + + AF 2 2 ' + , (27) 



where 



AF 2 ' e = lim — — [lndet N°' £ (n) + (n - 1) lndet JV 1,e (n)] , (28) 
AF 2 ' e = lim — I— y"(c£- C^" 1 ) lndet iV fc ' £ (n). 

For k > 2 one has (C* - C*" 1 ) ~ 2*=* so that 



1 + °° 9t — i 

^■'=m^ iri ^'" detMt ''' (29) 

There is a subtlety in the computation of AF± ,+ : as the limits of N 0,+ (n) and N 1,+ (n) when 
n — » are the same, one may be tempted to say that in this limit we have a unique family of 
eigenvalues of multiplicity re, and so AF 1 ' + = l/(2/3iV) In det M 1,+ . It is what || did, but it is 
wrong. Actually there is a factor 1/n to take into account, so that one also gets the contribution 
of the derivatives 



AF?<+ ' 



1 2/3N 



dlndetN°>+(n) dlndetJV 1 - 



n 



dn dn 



+ In det M 1,+ 

n=0 



(30) 



It turns out that the extra term is zero when (3 — > +oo, but it is not trivial (see appendix 
A.l ), We also show en passant in this appendix that detM 1 '" 1 " has a non zero finite limit when 
(3 —* +oo so that we eventually agree with || on the fact that 

AF 2,+ = 0. (31) 

As far as the computation of AF 2 ( |29| ) is concerned, the strategy of || consists into trans- 
lating the infinite dimensional matrices M k,e into more tractable integral operators. If (c p ) is 
an eigenvector of M k,£ , then 

/(*) = £ ^V^c q e( k+ ^- G ^ 2 (32) 



g=0 



gi 



is an eigenfunction, with the same eigenvalue, of the operator 

M k ' £ (x, y) = 5(x -y)- e(-l) k A k (x, y), (33) 
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A»(„ „, = 2»p ( -<M±m + ^ + „)) + f J-^' (34) 



where 

, /)=0 p!(2fc+p-l)! 
and reciprocally. 

The article [|| , on the basis of numerical discretization and diagonalization of these opera- 
tors, argues that the values of det M fc ' £ (T) plotted versus Tin A; fall onto two universal curves, 
depending on the parity of k: 

de t MM (T) .{™ I"-™ . (35) 

A more accurate numerical analysis showed us that this happens only when k — > +00, Tin A; 
being kept fixed: /+ and /_ are limit functions (see fig. [[]). Happily this does not change the 
conclusion that, in the limit (3 — ► +00, 

1 r+00 

AF "' £ = ~2N i dt[lnf£it) ~ ( 36 ) 

In ||, this integral is performed by fitting the numerical curves of /+ and /_ by smooth 
functions, which yields an estimate flawed by a rather rough uncertainty: 

AT 2 2 •+ = — (0.47 ±0.05). (37) 
As a consequence, AT itself is known with a bad precision: 

AT = —(—0.13 ± —0.05). (38) 

It is possible to improve on this. On can explicitly compute the limit of the operator A k (x, y) 
when k — > 00 under the restriction that t = T\nk remains fixed. In this case one also has 
j3 — > +00. Let us recall that in this limit 

G{1) = G(f3l) where G(x) = ln(l + e 2x ). (39) 

So we can write 

A k (x,y) ~ — ? _ ———= /(fc,exp(x + y)), (40) 

/3V(1 + exp(2s//?))(l + exp(2y//?)) 

where 

/(fc,z) = z k y - Liii - , (41) 

The eigenvalues of A fc (x,y) are the same as the ones of the operator 

Ink . k/ lnk Ink . 2 lnfc . ... , jlri . 

— A fc (— u,— v) = /(fe >e xp — (u + «))). (42) 

1 t t + exp(2u))(l + exp(2wj) c 

Let us define 

g(/c, w) = /(/c, exp(w In fc)). (43) 



In appendix A. 2 we show that when k — > +00, g(k, w) — > 0(w — 2) where is the Heaviside 



function. So the operator we have to diagonalize is 
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f "^* + + + + + * + + + + + + 



1.12 - sc" 
■ 

1.1 - 
1.08 -* 
1.06 - 



k=2 
k=3 
k=4 
k=5 
k=6 



k=2 
k=3 
k=4 
k=5 
k=6 



i □ * ■ — e^- 



Figure 1: The values of det M k,e (T) denned in (^) plotted versus Tlnk. Above, detM fc ' k 
odd and detM fc ' + k even. Below, detM fc, ~ k odd and detM fc,+ k even 



S 



8(^-2). (44) 



V(l + exp(2n))(l +exp(2u)) t 
It is the same as diagonalizing 

H t (u,v)= , 2 - @(u + v). (45) 

V ' y/(l + exp(2(u + t)))(l + exp(2( W + t))) 



We have /+(i) = lndet(7 — i7i) and /-(£) = lndet(7 + Ht). The correction pq ) reads 

9 J. 1 /' + 00 

Ai? 2 ' =7777/ d* [ln(det(7 + H t )) - ln(det(7 - 77 t ))] 
2iV Jo 



1 +2? /, 



±y^n±L (46 ) 

iV ^ 2p + 1 v ; 



where 



/•+oo /*+oo r 

I p = dtTrHf = / dt dui...du p Ht(ui,u 2 )...H t (up- 1 ,u p )Ht(u p ,u 1 ) (47) 

= 2° f °° dt[ dui ... dUp ( e^ + ^-.-e^-i + ^e^ + uO 

Jo 7 p (l + exp(2(t + -ui)))...(l + exp(2(t + -u p ))) 

(to derive (^) we used the identity lndet = Tr In and expanded ln(det(7 ± 77^)) in power series 
of fl*). 

Note that the operators Ht have positive and negative eigenvalues. The changes of variables 
Xi = exp(— 2ui) and C = exp(2i) yield 



-oo 



I p = dtTrHf = I ^/-^.../_J^E_ (48) 



f+ oc dC f dxi f dx p 
Jo 1 h 2C J xi + C" J Xp + C 1 

where the integration with respect to to be performed over the domain defined by 

Vi,Xi > and XiXi+i < 1, < 1. 

Unfortunately we were not able to compute analytically I p for a generic p. We succeeded 
in computing exactly the four first terms, and we got an estimate of the fifth one by numerical 
integration: 

C(2) 

h = ~ 0.411234, (49) 

h = ~ 0.601028, 

J 3 = ru 0.202936, 

16 

J 4 = 4C(5) - ~ 0.193102, 

7 5 ~ 0.137098. 

A truncated summation of (|4q) up to the third term gives the following lower bound 

AF 2 2 ' + > 0.506298/iV, (50) 
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Figure 2: (L min - tt 2 /12)N versus 1/N. The dashed line is the fit (pi) 



which is compatible with the estimate of Mezard and Parisi (|37|) . Adding AFi (|15|) , we get 

AF > -0.0947301/iV. (51) 
One can also try a Pade summation to compute AF 2 ' ■ One gets 

NAF ~ I x + — - — = -0.076. (52) 
1 1573 

To check the validity of ([51]) , we carried out numerical simulations similar to the ones in |Q] , 
but averaging over more samples and implementing a variance reduction trick [18, jl9|| . We used 
the values N = 35, 50, 60, 75, 100, 125 and 200, with a decreasing number of samples, from 
1200000 downto 300000. The results for L m i n are plotted in fig. [2|. A quadratic fit 

/- vr 2 \ b c 

[ L ^n-^N = a + - + — 2 (53) 
gives a = —0.0346 ± 0.0066, which is compatible with Q5"T|), 



3 Bipartite case 

Now we turn to the assignment problem. To make the comparison with the non bipartite case 
easier, we prefer using a slightly different convention for the partition function, which amounts 
to a rescaling of 0: we set 

Z = exp ^-y/3L matc /j in ^ (54) 

all possible matchings 

instead of (0) . However the reader must bear in mind that in this case L m i n is TWICE the free 
energy density F/N = -l/(/3N) In Z . 
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Moreover, as we want to compare our results with 



we must take into account both 



possible distributions of the distances. In appendix |B.l| we sketch the derivation of the averaged 
replicated partition function both in the case of the flat distribution (/i = 0) and in the case of 
the exponential distribution (/z = 1): 



N 



II dX a dY a - exp --S[X a , Y a ] 



N 



(55) 



1 



( £ %r£(X*j P + Y 2 u/3 ) ] exp ( xl + Y c 

z a>/3 9 aU (3 



where 



X 2 + Y 2 

S[X a , Y a ) = J2 — ~ - 2 ln z \X a ~ iY a ] - 2 lnz[X a + iY a \. 



9a 



(56) 



The thermodynamical limit of F/N does not depend on fj,, but its correction in 1/N does. 
One can look for a saddle-point of the particular form: 



X? = X s / {a) and Y°? = 0. 



X sp satisfy the equation 



X n 



2g Q 



9 In; 
~dX n 



(57) 



(58) 



which is exactly the same as the one for the non bipartite case ©. Hence 



X s / = Qf. (59) 

The free energy in the thermodynamical limit is the same as in the non bipartite case (|14|). 
Like in the non bipartite case, the 0(1/N) correction to the free energy contains the terms 
AF 1 coming from the last line in (|55|), and AF 2 coming from the Gaussian fluctuations. 
One has AF 1 = AF/ + /iAFf, with 



AF, 



2nNP^ 9 2 aU , [Q ' S 



and 



AF,} 



sp\2 



(60) 



(61) 



We see that AF/ is twice the one computed in the case of the simple matching problem ([T^) . 
AFg, computed in appendix B.2| , happens to be the opposite of AF/. So 



AF 1 = (/x - 1) 



03) 
N ' 



(62) 



which demonstrates our assertion (|5|), confirmed by numerical simulations [pTCj ] . 

As far as AF 2 is concerned, it is easy to generalize the computation of the non bipartite 
case. We have 



AF 



1 



2(3Nn 



In 



det det T ai 



(63) 
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where are the matrices of equation ( |T7| ) . 

n 2~l~ 2-f- 2 2 

Thus we can write AF Z = AF^' T +AF 2 '" r +AF 1 ' +AF 2 ' , where the different contributions 
are defined in We have Ai 7 ^' -1- + Ai 7 ^' - = (see eq. (|36|)) and we know that A.Fi' + = 

@. 

There is subtlety however in the computation of Ai 7 ^ ' , as pointed out in ||: we run i n t° a 
problem because T~p has some zero modes. This is actually no surprise: it arises from the fact 
that the action (pq) is left invariant under the transformation 



Xa + iY a (X a + iT„) exp i^9 fl , (64) 

\ aea J 

where 9\,...0 n are real angles. The zero modes are the n Goldstone modes of this invariance. 



As a consequence the saddle-point (|57|) is not unique: there is a n dimensional hypersurface of 
degenerated saddle-points parameterized as follows 



X? + iY? = X; p (a) exp (i J20a) < Bi < 2vr. 



(65) 



The kernel of T a/3 is spanned by the n vectors £j of components 



which have the replica % distinguished. So N k '~(n) for k > 2 has no zero mode. Only N 0, ~(ri) 
and A rl,_ (n) have a zero eigenvalue each. Thus Ai 7 ^' - is to be computed leaving aside the zero 
modes, and one has to take into account a new contribution Ai 7,3 to the free energy corresponding 
to the volume of the orbit ( |65|) . 

As far a the computation of AFi ' is concerned, we refer to appendix (IA.ll), where we 

2 + 2 — 

showed that Ai 7 ^ ' = 0. It is easy to see that the proof is exactly the same for Ai 7 ^ ' : there 
are only some sign reversals (in particular one has 1 + A/% instead of 1 — in the denominator of 
(|7q))) and the exclusion of the zero modes (corresponding to k = 1 in (|83|)) has no consequence 



because the key property ( [79[) holds for each eigenvalue of I(x,y). So Ai 7 ^' = 0, and 

AF 2 = 0. (67) 



Let us now turn to the computation of the volume of the hypersurface defined by (pB]). It 
is where Mezard and Parisi made a mistake: they computed this quantity without taking into 
account the fact that they carry out the diagonalization in another system of coordinates. To 
make things clearer let us rewrite (j55|) as 



Z n = Jl[ dU a dV a ^ exp (-^S[y/thP a , Vg^V a ]^j (68) 
x exp (-1 £ Ml(^ U7 + V^)\ exp (-^9^ a + V*)) . 

It is in the variables (U a , V a ) that Mezard and Parisi have chosen to diagonalize: indeed the 
matrix is half the Hessian matrix of S[y/g^U a , ^/g^V a \. So the volume of the hypersurface 
is to be computed in these same variables, not in (X a ,Y a ) as they did. The saddle-points 
coordinates are related by 
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The computation of the correct volume is done in appendix |B.3j . We find 

AF 3 



1 



(69) 



(70) 



Collecting the pieces we get the expected result (|) 



4 Conclusion 

By fixing the mistakes made by [|| in the computation of the 0(1/N) correction to the mean 
minimum length in the random assignment problem, we removed any inconsistency among the 
corpus of results about this problem. This gives further evidence that the replica approach, in 
its simplest symmetric ansatz, exactly solves the problem, and remains a valuable tool to gain 
insight on such quantities as this finite size correction which, despite recent dramatic progresses 
in the rigorous approach 0, still resist a mathematical treatment. 
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APPENDIX 



A Non bipartite case 

A.l Contribution of the derivatives in AFf ( |30| ) 



Our purpose here is to show that 



lim 

/3^+oo 2(5N 



dlndetiV°'+(n) din det N^ + (n) 



dn 



dn 



n=0 



(71) 



equals 0. 

We start from the equations ( pop and (p2[). Writing = 1 + ^ + o(n) it is easy to see that 
N 1,+ (n) — N 0,+ (n) = nA 1 + nA 2 + o(n), where Ai and A2 are the following infinite dimensional 
matrices 



-l)o +1 ( p + q -l)\Q p+q 
(jp -!)!<?! g p+q 



= ^ q V V.J ^Vftft (72) 
A 2 = (-1)' i +oc . (7;J) 

It is more convenient to use an integral operator formalism. We set 
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y ) = i ex P (- m+m^ K[x + y)t (74) 

where K is the function defined in ( |l3|). 

If (c p ) is an eigenvector of (|2l|) for the eigenvalue 1 — A then / is an eigenfunction of 

I(x,y) for the eigenvalue A, where 

m = E M^vw^-^)/ 2 . (75) 

One can easily check that the operators corresponding to the matrices A 1 and A 2 are re- 
spectively 



A\x,y) 
A 2 (x,y) 



e -G(x)/2 e G(y)/2 
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-G{x)e 



r+oo 

/ K(t + x)e- G{t) dt, 
Jy 



-G(y)/2-G(x)/2 



(76) 
(77) 



Thus ( |71[) can be evaluated by the standard result of first order perturbation theory 

r = _J_ V (WllMM, (78) 



2N ^ 0(1 - A fc ) 

where the \k) are the normalized eigenvectors of I(x,y) (|74|), and the A/j are the corresponding 
eigenvalues. 

Below we show that 

(fc|A 1 |fc) + (fc|A 2 |fc) = Q (79) 

for each k when — ► +oo, and that the A^ have finite limits, different from 1, so that T = 0. 
Let us consider an eigenfunction / of I(x,y): 



dyl(x,y)f(y) = Xf(x) 



We make the substitution 



f(x)e G ^ 2 = P(- 



1 + exp(2x//3) ; ' 

In the — ► +oo limit we can use (p{3) and after some changes of variables we see that 



(80) 



(81) 



can be restated as 

Vt> € [0, 1], 



1 — u 



0^((l-u)(l-v) 



XP(v) 



i.e. 



uv 

1- V 1—U 



(82) 



because K(j3z) = -1 if z > 0, otherwise. Note that P(0) = 0. 

The eigenproblem (82) is independent. We found that its eigenvalues are X k = (—l) k /k, 
k = 1,2,... + oo. The corresponding eigenfunctions are polynomial of degree k: 



P k (u) = ^2 a Pjk u p with a P)k 
p=0 



k\k>-i)...(k>-( P -iy) 

{ ' (p!)2 



(83) 
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For the computation of (k\A 1 \k), it is simpler not to use this explicit form. By derivation of 
), we get 



P(l-v) = XvP'(v) 
Combining (^), @ and (§0|), we have 



A 1 (x,y)f(x)dx = -Xe G ^ 2 / f(t)e~ G W 2 



+00 



so that 



A 1 {x,y)f(x)f(y)dxdy = A^ ^ du-^K C dv - P ^ 



This can be further simplified: 



4 Jo u(l - u) Ju 1 — v 



-a 2 — p(u) 



4 7o it(l — it) 



P(l-u) by (84). 



(84) 



(85) 



(86) 



1 , P(u)P(l-u) 

du r 

u(l - u) 







1 /ll 

duP(u)P(l -«)(- + 



u 1 — u 



so that 



2A f 1 duP{u)P'{u) 
Jo 

AP(1) 2 , 



A 1 (x,y)f(x)f(y)dxdy = -A 3 ^P(1) 2 . 



(87) 



The function / is a priori not normalized so that the above quantity is to be divided by 

r 1 . 



f(x) 2 dx = 2 J du Y 



u 



We eventually get 



(t|A 1 |t> = -A 3 |p(l) 2 



1 — u 



-1 -1 



Now we compute (k\ A 2 \k): thanks to (|T^), 

2 



| dt/(x,t)e- G ^/ 2 e- G W/ 2 



(89) 



(90) 



(91) 



hence 



A 2 (x,y)f(x)dx = A / dtf{t)e~ G ^/ 2 e~ G ^' 2 , 



(92) 
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Figure 3: The contour C in the complex plane 



& 2 (x,y)f{x)f(y)dxdy 



dtf(t)e~ G ^ 2 



(93) 



2 Jo 1-u 



i 2 



A 3 ^P(1) 2 



-i 



and (k\A 2 \k) = +A 3 §P(1) 2 £ du^- , which ends up the proof. 
A. 2 Limit of g(k,w) when k — > +oo 

We start from an integral representation of g(k,w) defined in (E3h : 



g(k,w) 



^kw In A; 



^\PgPt« In fc 



p=0 
„fcu> In 



2vr 



p!(p + 2A; - + k 

k ( dz [ +C °dxe s ^\ 
Jc Jo 



where 



S(k, z, x) 



k" 



-z - k(x + 21n(-z)) H e 



(94) 
(95) 



because (p + k)- 1 = J +o ° dx e^ p+k ^ and {p + 2k- l)!" 1 = i/(2n) f c dz e - z {-z)-^ +2k \ C is a 
contour in the complex plane such as illustrated in fig. 3. 
The stationarity equations read 



dz 

— k 

dx 

so that there is a movable saddle-point at 



k k % 



z 



0, 



(96) 



x sp = (w — 2) In k 



-sp 



-k 



(97) 
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A. 2.1 The case w > 2 

In this case the saddle-point (|97| ) is inside the range of integration. We compute the Hessian at 
this point: 



So, when k — ► +00, 
dxe s[k ' z ' x] ~ 



#5 

at 2 

cte 2 
dxdz 



dx exp 



e 5sp exp 



2— + 2— e - * = 0, 

z 

A: 11 



z 2 z 3 



(98) 



-A;, 



A; 1 



1. 



5 sp - - [x - (w - 2) In k] 2 + [x - (to - 2) In jfe] (21 + fc) 



2A- 



(2 + A;) 2 



Then we perform the integration with respect to z with z + k = —it: 



OO f 2 

/ / dee~zk 
00 



Given that S 



sp 



= —i\Fhxk. 
-kwlnk, we eventually get 

J J dxdze s{k ' z ' x) ~ -27rie- fcwlnfe . 



It follows that 



(99) 



(100) 



(101) 



g(k,w) 



*k— >+oo 



1. 



(102) 



A.2.2 The case w < 2 

In this case the saddle-point ( p7| ) is outside the range of integration. The integral is dominated 
by < x <C 1, where 



Jfe* 



A: 6 



S(k, z, x) = -z- 2k ln(-z) H (k H )z + 0(z 



so that 



exp 



-z-2fcln(-«) + V 



A: + 



(103) 



(104) 



We just have to look for the saddle-point of T,(k, z) = —z—2kln(—z)+k w /z. The stationarity 
of £ with respect to z reads 



2k k w 
1 s- = 0. 

z z z 



(105) 



There are two candidates saddle-point: 
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zf p = -k± Vk 2 - k w . (106) 

It is easy to see that, when k — ► +00, ) < whereas T,"(z~ ) > 0. So on a contour of 

the shape of A.2 the right saddle-point is z~ p . We have £(&, z^,) = — 2k In A; + 0(A;), so that 



/ dz / dze 5 **'*'*) = O e - 2fclnfc +oW . (107) 
Jc Jo L J 



/c JO 

Remembering (|94|), it follows that 

g(k, w)=0 r e fc(»-2)infc+o(fc)l _ ( 1Q8 ) 
So g(fc, to ) goes to when k — ► +00. 

B Bipartite case 

B.l Computation of the averaged replicated partition function 

We have two sets of N points each. We introduce the occupation numbers rnj = or 1, which 
are constrained by 

N N 

Viel,...N, ^riij = J2 n ji = 1 - ( 109 ) 

3=1 j=l 

The length of the matching associated to a choice of the is L({mj}) = J2i,j n>ijkj- The 
partition function (|34|) reads 



V 



X 



["Jexp li/xi I l-^jnj-ij j exp ( y >X "''A' ] > 



where we enforced ( 109 ) using an integral representation of the Kronecker symbol 5(p) 



fo dX/(2n)e i P x . It follows that 

z n = J d[\]d\fi] e Xp ^E A « a + ft a )n( 1 + E ex P (" f PW^i " i E W + M?)) J(iii) 

where d[A] is a shorthand notation for Yla=i Tli=i dX%/(2ir) (see the beginning of section || for 
any precision on the other notations). 

As we are interested in the subleading contribution to the free energy, it is important, when 
averaging over the disorder, to specify the distribution of the lengths we are considering: 



^ Np ^ /2 = ^9p-^4 + o( W 2) ( 112 ) 

where fj, = 1 in the case of the exponential distribution, in the case of the flat distribution. 
Thus 
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Z» = J d[\]d[n) exp (i E E X i + Mf) II (l + 2 ^ - 4 ^§l + °(^2 )) > ( 113 ) 



where 



a \ aGa / 

a \ aGa / 

Now we write 

11(1 + 2^-4^,) = exp(|E^)-p(-^E^) (US) 

x ex P " 4 t^ E ^ 

Plugged into (113), this gives 

W = j d[A]dMexp^EE A " + ^) ex p(^E^(4 + y')) (H6) 

X6XP (~2^2 E^f 7 (4u 7 +y a U 7 )) eX P -^E^^+^a) ' 
V a,7 / V a / 



where we set 



E ex P M E X i ) = ( x * + ^«)/ 2 ' 

i \ aGa / 

E ex P ( ~ { E Mi ) = i x a - iy a )/2- 

i \ aGa / 



(117) 



In ( 116[ ), YlijTjj gives the only contribution X)a, 7 because the other terms vanish when inte- 
grated (remember that for p integer, J d\e ipX = unless p = 0). 
Using well known properties of Gaussian integrals, we finally get 

W = [ d[X}d[fA exp (i EE X i + I*) I II dX a dY a ^- exp f-y E ^ ± ^) (118) 

X exp ( E ^^a + ^a?/a J exp ( - ~ E ^(^7 + ^U 7 ) J 6X P l-fi E *a + ^« j • 

\ a J \ Z a, 7 &*U 7 / \ a / 

Eventually, expressing x Q and y a as functions of the Af and /if one gets 

/jy / at y2 i y^2\ 

IJ dX a dF a - exp -- E z[X a - iY a ] N z[X a + iY a ] N (119) 

x exp f - \ E ^ + yJ U7 )) exp f- M £ j£ + *2 

V z q, j7 5aU 7 / V a 
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B.2 Computation of AF 2 X (|6lD 
Let us recall a result of [01 



2 r+oo e lp 



oo 



-G(Q 



Using ( p9| ) we see that, when n — > and /3 — > +oo, 
1 



E«2 



1 +°° c_np-i 



(120) 



(121) 



iV/3 2 



sp 



(_ lr l 



e p/ e -G(0 



2 f+°° 



dZG(Z) / (ite~ G(i) by (11) 



2 r+oo 

— du ln(l + «) ln(l + l/u) by 
iv Jo 



2N 

B.3 Computation of the volume of the hypersurface of saddle-points (|55|) 

The volume is J 27r d#i . . . J 27r y/detg where 



Safe 



E 



dd a de b + de a de b 



(122) 



<7 has a very simple structure: all diagonal elements are equal to go, all non diagonal elements 
to gi, with 



yp 



90 



P =i 



P =i 



gi = E(^ p ) 2 ^:2~^o/?E(- 1 )V(Q^) 2 + 5o. 

p=2 p=l 



(123) 
(124) 



Using fl59|) and (|120| ), one gets 

+°° O /-+oo „ip 



p=i 

" + 00 



oo (p-1)! 



-G(0 



(125) 



/ + 0O 
dlG'{l)e- G ^ = 2 
-oo 



and 



91 



+oo . 

go = 2j2(-V P PQp* 



+ 00 „lp 



dlG"(l)e- G ^ 



(p-1)! 
2 



(126) 



The computation of det g gives 
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so that 



det = [g (n) - 9l (n)]"" 1 ^^) + (n - l)ffi(n)], 
31 



\/det^= 1 + ^ 



In (50 - 3l] 



on 



(127) 



(128) 



_3o - 3i 

(here again one must be careful when deriving this result that there might a priori be some 
contribution of dg^/dn or dgi/dn). 

Thus the contribution to the free energy is 



AF 



1 



2N(3 
1 

"2N' 



31 



3o -3i 



+ ln(#o - 31 ) 



(129) 



Mezard and Parisi §] had expressions ( p3|) and (jl2|) with (jy , F p s P) instead of (C/ p s P, F p sp ), 
which made them find a wrong AF 3 = — 7r 2 /(24iV). 
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